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ABSTRACT 

We present a number of statistical tools that show promise for obtaining information on turbulence 
in molecular clouds and diffuse interstellar medium. For our tests we make use of three-dimensional 
512^ compressible MHD isothermal simulations performed for different sonic, i.e. A4s = Vl/Vs^ where 
Vl is the injection velocity, Vs is the sound velocity, and Alfvenic Ma = Vl/Va^ where Va is the 
Alfven velocity, Mach numbers. We introduce the bispectrum, a new tool for statistical studies of 
the interstellar medium which, unlike an ordinary power spectrum of turbulence, preserves the phase 
information of the stochastic field. We show that the bispectra of the 3D stochastic density field and 
of column densities, available from observations, are similar. This opens good prospects for studies 
of molecular clouds and diffuse media with the new tool. We use the bispectrum technique to define 
the role of non-linear wave- wave interactions in the turbulent energy cascade. We also obtained 
the bispectrum function for density and column densities with varying magnetic field strength. As 
expected, a strong correlation is obtained for wave modes ki = /c2 for all models. Larger values of 
Ais result in increased correlations for modes with ki ^ k^. This effect becomes more evident with 
increasing magnetic field intensity. We believe that the different MHD wave modes, e.g. Alfven and 
magneto-acoustic, which arise in strongly magnetized turbulence, may be responsible for the increased 
correlations compared to purely hydrodynamical perturbations. In addition to the bispectrum, we 
calculated the 3rd and 4th statistical moments of density and column density, namely, skewness 
and kurtosis, respectively. We found a strong dependence of skewness and kurtosis with M^s- In 
particular, as Alg increases, so does the Gaussian asymmetry of the density distribution. We also 
studied the correlations of 2D column density with dispersion of velocities and magnetic field, as well 
as the correlations of 3D density with magentic and kinetic energy and Jv[a for comparison. Our 
results show that column density is linearly correlated with magnetic field for high M-s- This trend 
is independent of the turbulent kinetic energy and can be used to characterize inhomogeneities of 
physical properties in low density clumps in the ISM. 
Subject headings: ISM: structure — MHD — turbulence 
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1. INTRODUCTION 

The interstellar medium (ISM) is known to be highly 
turbulent but the role of magnetic fields in the dis- 
tribution of turbulent eddies, as well as in the energy 
cascade, is still not completely understood. Numeri- 
cal simulations have proved to be important tools to 
study turbulent processes in th e magnetized ISM (see 
iBallesteros-Pared es. et all '2007). Stars are known to 
be formed on the denser regions of the ISM. However, 
the formation and survival of these dense cores is not 
completely understood, mostly because of the complex 
relationship between turb ulence, mag:netic fields, and 
self-gravity (see reviews by Lazarian & Cho 2004; Evans 
[L999; Elmegreen & Scalo 2004; Mckee & Ostriker 2007). 
In a rather simplistic view, highly supersonic turbulent 
motions generate shocks that evolve into dense struc- 
tures, which may be dispersed by turbulence. 

This picture may be very different if the magnetic 
field is taken into account. Observationally, line pro- 

^ Astronomy Department, University of Wisconsin, Madison, 
475 N. Charter St., WI 53711, USA 

^ Department of Physics Sz Astronomy, University of Louisville, 
Louisville, KY 40202, USA 

^ Astronomical Observatory, Jagiellonian University, ul. Orla 
171, 30-244 Krakow, Poland 

^ Nucleo de Astrofisica Teorica, Universidade Cruzeiro do Sul - 
Rua Galvao Bueno 868, CEP 01506-000, Sao Paulo, Brazil 



files indeed indicate ISM turbulence is supersonic, but 
the determination o f the mag;netic field intensity is much 
more involved (see iFalceta-Goncalves et al.l l2003f ) . In 
spite of the extensive use of current techniques, e.g. po- 
lar imetry and Zeeman splitting measurements, the ac- 
tual ratio of the magnetic and tur bulent kinetic ener- 
g;ies is still a matter of debate (see iPadoan fc Norlundl 
120021 : iGirart. Rao fc Marrond 120061 ). Magnetic fields 
may play a key role in the star formatio n proc ess (see 
iMcKee fc Taj] 120021 : iMac Low fc KlessenI [2QQ3 . for re- 
view), e.g. providing clo uds with extra suppo rt against 
gravitational collapse (see lOstriker el al.]ll999l ). Strongly 
magnetized turbulence presents weaker shocks and de- 
creased densities compared to purely hydrodynamic tur- 
bulence. 

Because of the difficulties mentioned above, numer- 
ical simulations represent a unique method to un- 
derstand the evolution of turbulence and the role of 
the magnetic field in the creation of d e nsity struc- 
tures in the ISM. Vazquez- Semade ni et al.l (|l997l ) stud- 
ied the Larson-type correlations of density and mass 
with cloud radii, using two-dimensional MHD simu- 
lations. It was shown that the topology of density 
is filamentary and its probability distribution function 
(P DF) is, in most cases, log-normal. Further stud- 
ies (IScalo et all[l99 8: Vazquez- Semadeni fc Garcia 20^; 
lOstriker el al.ll2QQlii) confirmed these results (see Beres- 
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nyak, Lazarian & Cho 2005). Statistical descriptions 
of turbulence are valuable as they constrain the physi- 
cal properties of the system. Several techniques involv- 
ing numerics have been p roposed f or velocity studies 
(see Lazarian & Pogosvan 2000; Lazarian fc Choi l2004l : 
[Lazarian fc Pogosy an 2 006 , 2 008 . ) and have be e n suc - 
cessfully tested in works such as iLazarian et al.l (|200ll ): 
lEsauivel et aH (|2003D : iPadoan et al.l (|2006[ ). These stud- 
ies are aimed in providing the velocity spectrum of the 
turbulence. It is important to note that some of the 
earlier results obtained for velocity spectra studies using 
velocity centroids, are likely to be in error for high Mach 
number turbulence in molecular clouds in view of the re- 
cent insight into the properties of velocity centroids (see 
Esquivel & Lazarian 2005, Esquivel et al. 2007). 

In the past, density studies have also provided insight 
on astrophysical turbulence. Although density is a less 
direct measure of turbulence compared to velocity, it is 
readily available from the observations of column densi- 
ties. The index of the density spectrum gets shallow with 
the sonic Mach number Ais = Vl/Vs^ where Vl and Vs 
are the injection and sound velocities, respectively, as 
shown in MHD simulations by Beresnyak, Lazarian & 
Cho (2005). This was also confirmed to be also true for 
hydrodynamic simulations in the follow-up simulations 
in Rye & Kim (2006)^. The shallow spectra of density, 
in fact, corresponds to observations (see Lazarian 2008 
for a review) of shallow density spectrum of column den- 
sities measured for molecular cloud studies. 

The complexity of turbulence calls for the simultane- 
ous use of different statistical measures. Recently, Kowal, 
Lazarian, & Beresnyak (2007) (henceforth referred to as 
KLB) tested how different measures, including PDFs, 
power spectrum, skewness, kurtosis, and higher order 
She-Leveque exponents vary with and Ma using an 
extensive set of MHD simulations. They also provided a 
study of the variations of topology, using genus, ^ as well 
as anisotropics of density structures, using correlation 
functions, with Aig and Ma- In KLB, emphasis was 
made on comparing the results available through syn- 
thetic observations of column density and the underlying 
properties of 3D density. This made this study easily ap- 
plicable to observations, e.g. to Wisconsin H a Mapper 
(WHAM) data (see Hih et al. 2008). 

Despite the number of past works, the potential of den- 
sity studies is far from being exhausted. For instance, 
when the energy cascade process and the distribution 
of density structures are analyzed from the power spec- 
tra of two-point correlations, phase information is lost. 
On the other hand the bispectrum, or three-point statis- 
tic, measures the magnitude and the phase of the cor- 
relation of signals in Fourier space. As a consequence 
of being a more informative measure, it can be used to 
search for nonlinear wave-wave interactions and charac- 
terize the turbulent regime. The bispectrum has been 
wide ly used in cosmology and g:ravitationa l wave studies 
(|£ry||l998' 'Scoccimarro''2000';'Liguori et al.""2006) and for 
the characteriz ation of wave-wave interactions in labo - 
ratory plasmas (jlntrator et al.lll989l : iTvnan et al.ll200ll ). 

^ The first mentioning of the effect can be traced back to Padoan 
et al. (2004) study. 

^ Genus was earher briefiy discussed in the context of interstellar 
medium in Lazarian (1999), Lazarian, Esquivel & Pogosyan (2001) 
and Lazarian (2006). 



and was suggested in iLazarianl (1 19991) to be applied to 
ISM turbulence (see also ILazarian. Kowal fc BeresnyakI 
I2008D . In this paper we apply the bispectrum to density 
and column density maps for different turbulent regimes 
in order to test the usefullness of the tool for interstellar 
studies. 

In addition, this work extends the study in KLB in sev- 
eral respects including the use of higher resolution simu- 
lations, further investigation of the dependence of Mach 
number on the skewness and kurtosis of densities, and 
studies of 3D correlations, which were not attempted in 
KLB. For instance, we study the correlations between the 
physical parameters (velocity dispersion, magnetic field 
intensity, etc.) with density and column density, and in- 
vestigate the correlation between velocity dispersion and 
LOS velocity, in order to explore their evolution within 
dense structures and their correlational dependence with 
turbulent regimes (namely the sonic and Alfvenic Mach 
numbers). We also provided the first interstellar medium 
related study of the bispectrum. In particular, we study 
the bispectrum of density and column density for sev- 
eral models of MHD turbulence. In § [2] we describe the 
numerical models of compressible MHD turbulence. In 
§ [3] we discuss the skewness and kurtosis of density and 
maps of column density. In § |4] we examine the corre- 
lations of Mach numbers, magnetic and kinetic energies 
with density and magnetic field and velocity with column 
density. In § [5] we present an analysis of the bispectrum 
of compressible MHD turbulence. In § 6 we discuss our 
results, followed by our conclusions and the summary. 

2. NUMERICAL SIMULATIONS 

We used a second-order-accur ate hybrid essentially 
nonoscillatory (ENO) scheme (see lCho fc Lazarianll2002T ) 
similar to that shown in KLB to solve the ideal MHD 
equations in a periodic box. 



dp 
dt 



V-(pv)=0, (1) 
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dB 

~dt 



V X (v X B) = 0, (3) 



with zero-divergence condition V-B = 0, and an isother- 
mal equation of state p = c^p, where p is density, v is 
velocity, B is magnetic field, p is the gas pressure, and 
Cs is the isothermal speed of sound. On the right-hand 
side, the source term f is a random large-scale driving 
force (in fact, we drive only the velocity field). The RMS 
velocity SV is maintained to be approximately unity, so 
that V can be viewed as the velocity measured in units of 

1/2 

the RMS velocity of the system and B/ (47rp) ^ as the 
Alfven velocity in the same units. The time t is in units 
of the large eddy turnover time (~ L/5V) and the length 
in units of L, the scale of the energy injection. The mag- 
netic field consists of the uniform background field and 
a fluctuating field: B = Bext + b. Initially b = 0. 

We drive turbulence solenoidally, at wave scale k equal 
to about 2.5 (2.5 times smaller than the size of the 
box). This scale defines the injection scale in our mod- 
els in Fourier space to minimize the influence of the 
forcing on the generation of density structures. Den- 
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sity fluctuations are generated later on by the inter- 
action of MHD waves. Density structures in turbu- 
lence can__be_associated with the slow and fast modes 
fsee iLithwick fc G oldreich 2001; Cho fc LazarianI l2003l : 
iBeresnvak. Lazarian fc Chd 120051 ). We use units in 

1 /2 

which Va = ^ext/ (47rp) =1. The values of 5ext have 
been chosen to be similar to those observed in the ISM 
turbulence. The average RMS velocity 5V in a statisti- 
cally stationary state is around 1. 

We do not set the viscosity and diffusion explicitly 
in our models. The scale at which the dissipation 
starts to act is defined by the numerical diffusivity of 
the scheme. The ENO-type scheme s are c onsidered to 
be relatively low diffusion o nes (see iLiu fc Osher 1998; 
iLevv. Puppo fc Russo 1 1 19991 . e.g.). The numerical diffu- 
sion depends not only on the adopted numerical scheme 
but also on the "smoothness" of the solution, so it 
changes locally in the system. In addition, it is also a 
time- varying quantity. All these problems make its esti- 
mation very difficult and incomparable between different 
applications. However, the dissipation scales can be es- 
timated approximately from the velocity spectra. In the 
case of our models we estimated the dissipation scale kj^ 
at 30 for high resolution. 

We present 3D numerical experiments of compress- 
ible (MHD) turbulence for a range of Mach numbers. 
We understand the Mach number to be defined as the 
mean value of the ratio of the absolute value of the lo- 
cal velocity V to the local value of the characteristic 
speed Cs or Va (for the sonic and Alfvenic Mach num- 
ber, respectively). We divided our models into two groups 
corresponding to sub- Alfvenic (^ext = 1-0) and super- 
Alfvenic (^ext = 0-1) turbulence. For each group we 
computed several models with different values of pressure 
(see Table [T]). We ran 6 compressible MHD turbulent 
models, with 512^ resolution, for t ~ 5 crossing times, to 
guarantee full development of energy cascade. We show 
in Figure [1] the full development of global quantities mag- 
netic, kinetic, and total energy. This figure shows that 
1 dynamical time unit (time normalized by the crossing 
time) is sufficient for saturation of energies. Since the 
saturation level is similar for all models and we solve the 
isothermal MHD equations, the sonic Mach number is 
fully determined by the value of isothermal sound speed, 
which is our control parameter. The different initial con- 
ditions led to three values of sonic Mach number, ~ 0.7, 
~ 2, and ~ 7, and two Alfvenic Mach numbers, ~ 0.7 
and ~ 2. The models are listed and described in Table [H 

Ultimately, we would like our numerical results to 
be directly comparable to observations. Thus, an ob- 
server might be curious about what density, tempera- 
tures, units, and scalings are appropriate for these mod- 
els. These are isothermal scale- free simulations. For the 
situations when gas can be assumed isothermal, they can 
easily be rescaled for any parameters (i.e. Temperatures, 
scales, and densities, see Table [2]) of the media con- 
cerned. In Table [2] the abbreviations stand for the fol- 
lowing six idealized interstellar environments: cold neu- 
tral medium (CNM), warm neutral medium (WNM), 
warm ionized medium (WIM), molecular cloud (MC), 
dark cloud (DC), reflection nebula (RN), and photodis- 
sociation region (PDR). We list suggested values for uh 
(Hydrogen density), T(k) (gas temperature), x (starlight 



intensity relative to the average starlight background), 
and other properties of these phases. These properties 
should be scaled in the models as dictated by the statis- 
tics used. An example of these scalings can be found 
in Hill et al. (2008), who used scaling relationships to 
set temperature, length, and density of simulated data 
to match the WIM. They explored PDFs for the WIM 
using data taken from the WHAM and synthetic data 
using these simulations. They found that scaled models 
with certain Mach numbers matched the WIM very well. 
Similar examples of simulation scaling can also be found 
in Falceta-Goncalves et al.(2008). In general, these simu- 
lations should be scaled accordingly to the characteristics 
of the media studied. These characterist ics are listed in 
Table [2] as shown in Drain e fc LazarianI (|199 8) . As for 
units, these are defined in terms of the speed of sound, 
which varies for a given Mach number. For A^5=7.0, 
C5=0.1, for Ms=2.0, C5=0.333 and for Ms=0.7, 0^=1.0. 

For reference, we include a cube visualization in Fig- 
ure [2] showing the X, y, z axes, the line of sight (LOS), 
and the direction of the magnetic field. We will refer to 
these directions throughout the paper. For example, "x 
column density" or "column density in the x direction" 
refers to the density cube being integrated along the x 
direction to make the 2 dimensional field. This descrip- 
tion is similar for the y and z directions. In the case 
of Figure [21 the line of sight is along the z axis and the 
magnetic field is oriented along the x axis. 

3. STATISTICAL MOMENTS: SKEWNESS AND KURTOSIS 

Skewness and kurtosis are defined by the third and 
fourth-order statistical moment. Skewness is defined as: 

If a distribution is Gaussian, the skewness is zero. Neg- 
ative skewness indicates the data is skewed in the left 
direction (the tail is extended to the left) while positive 
values imply that the distribution is skewed in the right 
direction (the tail is extended to the right). 

Kurtosis is a measure of whether a quantity has a dis- 
tribution that is peaked or flattened compared to a nor- 
mal Gaussian distribution. Kurtosis is defined in a sim- 
ilar manner to skewness, only is derived from the forth 
order statistical moment. If a data set has positive kurto- 
sis then it will have a distinct sharp peak near the mean 
and have elongated tails. If a data set has negative kur- 
tosis then it will be flat at the mean. Kurtosis is defined 
as: 




Skewness and kurtosis are statistical measures that 
can be used to characterize the properties of turbulence. 
They have been recently been used to characterize com- 
pressible MHD turbulence for simulations and synthetic 
observations (see KLB). In this work we test their results 
applying these calculations to cubes of finer resolution 
(512^) and further investigate how applicable this mea- 
sure is to observations. The interest on these measure- 
ments comes from that fact that the skewness and kurto- 
sis of column density and 3D density depends strongly on 
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Fig. 1. — Energy density vs. time, showing the evolution of magnetic, kinetic, and total energy for a A4a = 2.0, A4s = 0.7 model. Time 
is in units of the "crossing time," which is the time it takes the sound speed to travel through the cloud. It is defined as tcrossinq = — - — • 
Here we show that 1 time unit guarantees full development of the turbulent energy cascade. The difference in saturation energies is also 
seen in incompressible turbulent dynamo simulations, where it is of the order 30% (Cho et al. 2009) 



TABLE 1 

Description of the simulations - MHD, 512^ 



Model 


P 


^ext 


Ms"" 


Ma"" 


Description 


1 


1.00 


1.00 


0.7 


0.7 


subsonic &; sub-Alfvenic 


2 


0.10 


1.00 


2.0 


0.7 


supersonic &; sub-Alfvenic 


3 


0.01 


1.00 


7.0 


0.7 


supersonic & sub-Alfvenic 


4 


1.00 


0.10 


0.7 


2.0 


subsonic & super-Alfvenic 


5 


0.10 


0.10 


2.0 


2.0 


supersonic &; super-Alfvenic 


6 


0.01 


0.10 


7.0 


2.0 


supersonic &; super-Alfvenic 



• Sonic Mach number is defined as JVlg = (|v|/a}, where a is the sound speed and the averaging is taken over whole volume 

• Alfvenic Mach number is defined diS A4a = H^I/'^a), where va = \'^\/\/p is the local Alfven speed and the averaging is taken over 
whole volume 



TABLE 2 

Idealized Environments for Interstellar Matter 



Parameter 


DC 


MC 


CNM 


WNM 


WIM 


RN 


PDR 




10^ 


300. 


30 


0.4 


0.1 


10^ 


10^ 


T(K) 


10. 


20. 


100. 


6000. 


8000. 


100 


300 


X 


10^ 


0.01 


1.0 


1.0 


1.0 


1000 


3000 










0.0012 


0.1 


0.99 


0.001 


0.0001 


XM =n (M+) /nn 


10-6 


0.0001 


0.0003 


0.0003 


0.001 


0.0002 


0.0002 


y= 2n(H2)/nH 


0.999 


0.99 











0.5 


0.5 



From Draine &; Lazarian 1998. 



Ais (see KLB). Physically, this implies that one can de- 
termine M-s from examining the skewness and/or kurto- 
sis of the observational distributions. Also, by comparing 
the observed column densities to the synthetic ones and 
determining the turbulent regime of a given cloud, it is 
possible to understand the three-dimensional statistics. 

Statistics are obtained by applying Equations (|4]) 
and ([5]) to data cubes of density and column density. In 
Table [S] we show the values for the skewness and kurtosis 
of density for models with resolution of 512^, as well as 
for the synthetic observations, i.e. the column density 
maps (512^), for all models. We applied the moment 
statistics to 3 different snapshots for each model. Ulti- 
mately, we averaged together statistics from all 3 snap- 



shots to obtain the presented results with error bars de- 
termined using the standard deviation of the results for 
each model. These error bars are not the error on the cal- 
culated skewness of the models. Instead they represent 
the spread of skewness between the snapshots of each 
model, i.e. the standard deviation. From this data we 
can see how the skewness and kurtosis of these quantities 
depend on sonic and Alfven Mach numbers. 

The skewness of density and column densities in the x 
and y directions are strictly positive. We see the same 
trend in both density and column density in that the 
asymmetry grows with increasing Ais for both quanti- 
ties. As expected, small Ais niodels are more Gaussian 
(7^0), compared to supersonic models. The presence 
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Fig. 2. — Throughout the paper we make mention of directions x, y, and z as well as the line of sight (LOS). This figure shows an example 
of a subsonic, sub-Alfvenic 512^ data cube demonstrating one possible position of the line of sight (LOS) and mean magnetic field along 
the X direction in respect to the x, y, z axes. Here the LOS is orientated along the z axis, although throughout the paper we will investigate 
statistics for LOS along all three axes. 



TABLE 3 

Skewness and Kurtosis of Density and Column Densities Integrated along X, Y, and Z-direction 



Ms 


Ma 


P 












Skew 


Kurtosis 


Skew 


Kurtosis 


Skew 


Kurtosis 


Skew 


Kurtosis 


0.7 


0.7 




0.6±u-^ 


0_47±u.ui 






02±u.6 


-0.03±^-^^ 


o.oi±^-^ 


0.7 


2.0 


L0±0-2 


2.0±0-7 


Q37±o.oi 


_0.4±o.4 


0.18±0-04 


_0.7±o.6 


-0.06±002 


_0.2±o.2 


2.0 


0.7 


3.8±0-i 


20^3 


0,7±o.2 


0.8=^0-4 




2.i±o.4 




2±i 


2.0 


2.0 


4,0±o.i 


30±3 


L0±0-2 


l,4±o.4 


0.9±oi 


1 5=^0-^ 


0,6±o.3 


0.4±i 


7.0 


0.7 




200^30 




2±2 




3.6±0-9 




4±i 


7.0 


2.0 


10.1=^0-4 


300±30 


L6±0-3 


5±2 


L2±0-i 


2,4±o.9 


L2±0-i 


2±i 



of shocks, which result in high density structures, is en- 
hanced as A^s increases and these shocks in supersonic 
turbulence are the cause for the asymmetry in the den- 
sity and column density distributions. We also analyze 
the skewness of averaged data cubes as shown in Ta- 
ble m We average together points in a cube in order to 
smooth out far outlying data that might be affecting the 
overall skewness. From Table IH we can see that even af- 
ter averaging, the asymmetry in the density and column 
density distributions still increases with increasing Ms^ 
although not as sharply as the non-averaged case, due to 
the smoothing of outliers. 

The kurtosis of both density and column density is 
shown to be also higher for supersonic models, which 
implies sharper distributions. It is a consequence of the 
transfer of mass from the average values to the right tail, 
mostly due to shocks. Again, subsonic models present 
kurtosis around zero, i.e. similar to Gaussian, in agree- 
ment with KLB. Highly supersonic column density mod- 
els are systematically more peaked then subsonic models, 
following the trend noted from the density distributions. 

From this information we can see that the asymme- 
try of the distribution is greater for cases of high A^^, 



where shocks produce high density clumps. Just as in 
KLB, skewness and kurtosis both seem to be relatively 
unaffected by the strength of the magnetic field. In this 
way, we can determine by examining the skewness of 
observational column densities. For kurtosis, we see that 
subsonic models present nearly perfect Gaussian distri- 
butions for both densities and column densities. This 
implies that if density is perturbed weakly, it remains 
nearly normal. Both measures are strongly dependent 
on the sonic Mach number, making them useful statis- 
tics for observers to characterize interstellar turbulence. 



4. CORRELATIONS 



4.1. 



Correlation of three-dimensional fields 

It is generally accepted that isothermal simulations, 
such as the one presented in this paper, are reasonable 
models of molecular clouds, while non-isothermal sim- 
ulations better represent the many temperatures of the 
multiphase ISM. In order to study relationships between 
global parameters in molecular clouds or ionized gas, cor- 
relations are useful. At the most basic level, correlations 
can be used to determine a possible relationship between 
two quantities. Correlation plots can provide information 
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TABLE 4 

Skewness of Density and Column Densities Integrated along X, Y, and Z-direction averaged with a Gaussian beam of 100 

PIXELS 



Ms 


Ma 


7 Av. p 


7 Av. Ex 


7 Av. Ey 


7 Av. Ez 


0.7 


0.7 


0.7±0-i 


0.05±002 


0.3±0-2 


-0.04±0-2 


0.7 


2.0 




0.02±0-2 


0,2±o.3 


-0.1±0-2 


2.0 


0.7 


3.1±0-i 


0.63±oo4 


1.03±002 


0.9^0-2 


2.0 


2.0 


3.1±0-2 


0.8±0-4 


0.9^0-2 


0.5±0-i 


7.0 


0.7 


5.65^0-04 


Q_g7±0.02 


1.2±0-2 




7.0 


2.0 


5_Q±0.3 


1.5±0-2 


l.l±o-i 


l.OOi 



on the dynamical importance of kinetic and magnetic en- 
ergy in the evolution of density clumps in star forming 
clouds. In this section we study these correlations, pre- 
senting plots of magnetic energy vs. density, and Ma 
vs. normalized density. We also discuss the correlation 
of specific kinetic energy vs. density. In order to make 
comparisons to observations, plots of magnetic field vs. 
column density and velocity dispersion vs. column den- 
sity in directions parallel and perpendicular to the mean 
magnetic field are included. We also study the correla- 
tion between dispersion velocity for constant density vs. 
non constant density in order to see how density fiuctu- 
ations affect this observable quantity in our simulations. 
Correlations were preformed for the final snapshots of all 
models. 

The plots of magnetic energy vs. log density are shown 
in Figure [31 For supersonic cases we see that the mag- 
netic energy is correlated with density, namely, denser 
regions contain stronger magnetic fields, which is due to 
the compressibility of the gas and magnetic fields be- 
ing frozen in. This causes the magnetic field to follow 
the fiow of plasma if its tension is negligible. The com- 
pressed regions are dense enough to distort the magnetic 
field lines, enhance the magnetic field intensity, and effec- 
tively trap the magnetic energy due to the frozen-in con- 
dition. The magnetic energy in the super- Alfvenic model 
also reaches peaks higher than those in sub- Alfvenic cases 
due to a larger magnetic pressure in the latter case. For 
subsonic cases we see a very different trend as a result 
of incompressible turbulent fiows. Due to lack of shocks 
trapping the magnetic field in the density clumps, the 
magnetic energy is anti-correlated to density for sub- 
Alfvenic turbulence, which shows a narrow PDF peaked 
at the mean density. For super Alfvenic cases, the dis- 
tribution is more spread and no relationship is obtained. 

We examined the correlation of squared velocity (i.e. 
the specific kinetic energy {8k = V^) vs. the logarithm of 
density, but do not include plots since no clear relation- 
ship is seen. We observe a decrease in velocity for higher 
density, however here we are also interested in how den- 
sity and the kinetic energy density (i.e. = p ^ V'^ 
vs. p) relate to one another. For specific kinetic energy 
Ek vs. p, the distribution tends to isotropically cluster 
around the mean value of density. Therefore, while the 
specific kinetic energy might decrease with larger densi- 
ties, the kinetic energy density increases in a near lin- 
ear relationship for supersonic models. Subsonic models 
have minimum kinetic energies at values of log-density 
that vary less from the mean density then do supersonic 
models (-0.6 to 0.4 for subsonic and -2 to 2 for super- 
sonic). Compressibility in the subsonic models plays less 
of a role than in supersonic ones, and thus they reach 



their minimum kinetic energy with clumps that are less 
dense then supersonic models. Because density clumps 
are not prevalent in subsonic models, we see no clear re- 
lationship between density and kinetic energy, with the 
highest values for being around the mean of density 
and then falling off. 

Apart from the global Alfven Mach number, which is 
defined SiS Ma = ^ where Vl is the velocity fiow ^, 
which is computed by taking the magnitude of velocity 
Vl = \J ~^ ~^ '^z ^ (the Alfven speed), one 

can also determine the local Alfven Mach number, which 
is what observers can potentially measure ^. Ma de- 
pends on both the mean magnetic field and density p, 
which can result in differences between the global Ma 
of a cloud and a density clump within the cloud. In 
order to explore the relationship between density and 
Ma we present here the correlation of the log Alfvenic 
Mach number versus the log of , shown in Figure [H 
Although we use a log-log scale, the actual numbers of 
Ma are used for ease of reading. We obtained a general 
trend of Ma slightly increasing with density, especially 
for sub Alfvenic cases. This is because Ma has depen- 
dence on density: Ma — ^ — . Cases with 
globally sub- Alfvenic turbulence show local Ma contin- 
ually increasing with density until local values become 
super-Alfvenic. Cases that are globally super-Alfvenic 
reach local values that are higher then Ma = 2.0 which 
peak at the mean density. However, values of density 
larger then the mean yield decreasing values of Ma for 
the globally super-Alfvenic models, although they con- 
tinue to stay locally super-Alfvenic. In essence, we find 
that globally sub- Alfvenic models become locally super- 
Alfvenic with increasing density while the globally super- 
Alfvenic models remain so for all values of density stud- 
ied. This can best be explained by Mas dependence 
on B and the relationship between B and p. From Fig- 
ure [3l it is noticeable that the magnetic energy does not 
grow as effectively with density in sub- Alfvenic cases as 
it does in super-Alfvenic ones. Therefore, even in glob- 
ally sub- Alfvenic models, the very dense regions may be 
super-Alfvenic, in spite of the global regime of the cloud, 
since density clumps do not effectively trap the higher 
magnetic field. However, most density clumps lie around 
Ma = 0.7 and only the highest density clumps are able 

^ In this section, we denote global Alfven Mach number with a 
calligraphic A4a, versus the local Alfven Mach number which is 
denoted by Ma 

^ The estimate of the local Mach number is possible, for in- 
stance, using the Bz component of magnetic field available through 
Zeeman splitting and the velocity dispersion of clouds available 
through spectroscopic measurements. 
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Fig. 3. — The 3D correlation of magnetic energy vs. the logarithm of normalized density. The first row consists of super- Alfvenic cases 
while the bottom row is sub-Alfvenic. Images are ordered from left to right as supersonic to subsonic. Blue contours indicate regions of 
high data counts while red and yellow have lower counts. 



to become super-Alfvenic. Similarly for 7Uyi=2.0, past 
the mean density the Alfvenic Mach number begins to 
decrease. This is because the magnetic field increases 
as the density of clumps increase, causing the clumps to 
become more magnetically dominated and thus causing 
Ma to decline. As the models progress from subsonic to 
supersonic, we see that these effect becomes more pro- 
nounced. 

4.2. Synthetic Observational Correlations 

In order to make our studies applicable to observa- 
tional analysis, we provide correlational studies using col- 
umn densities. We performed this analysis for magnetic 
field and velocity dispersion in order to gain a better 
understanding of how these quantities affect turbulent 
gas. Interest in these measures is due to the fact that 
the magnetic field intensity parallel to the line of sight 
(LOS) is directly compared to Zeeman measurements, 
and perpendicular comp onents may be estimated using 
polarimetric maps (see iFalceta-Goncalves et al.l [2008). 
Regarding the velocity field, we focused on correlating 
column density to dispersion velocity along the LOS as 
obtained from line profiles. We also compare in Fig- 
ure [To] the correlation of theoretical velocity dispersion 
for constant density with actual velocity dispersion avail- 
able from synthetic observations. The goal of this is to 
study how density fluctuations affects the observables. In 
Figure [21 we present a cubic depiction of what is meant 
by the directions x, y, z and the LOS. We can place an 
observer at different points along the plane of the cube 
to look vertically through with parallel sightlines. This 
visualization utilizes every cell in the simulated cube. 
However, in reality we only have one direction of view- 
ing molecular clouds. The directions presented here are 
all measurable and mearly used to get a variety of col- 
umn density samples along different lines of integration. 



However, a cloud with a very specific viewing geometry 
might require integration along a sightline above or be- 
low a certain latitude (see Hill et al. 2008) and we do not 
specifically consider all these cases here. 

In Figure [5] we show the correlation of column density 
along the x-direction and the magnetic field component 
parallel to the LOS, i.e. parallel to 5ext- This correla- 
tion may be compared, for instance, directly to Zeeman 
splitting measurements. Similarly to the 3D correlation 
of Figure [3l column density and magnetic field along the 
LOS increase together for supersonic models due to den- 
sity clumps trapping the magnetic field. As expected, 
sub-Alfvenic models present steeper correlations due to 
higher mean magnetic field intensity. Correlations are 
near linear for column density values around the mean 
density for supersonic sub-Alfvenic models. This is due 
to the stronger field becoming entangled with high den- 
sity clumps. In super Alfvenic cases, the magnetic field 
is not as strong and hence density clumps will increase, 
yet the magnetic energy available for clumps to trap will 
level off. Subsonic models showed no correlation due to 
incompressible densities. 

A behavior similar to Figure [5] was obtained for the 
correlation of column density along x-direction and the 
magnetic field component perpendicular to the line of 
sight^, i.e. LOS ^- ^ext, shown in Figure [6l The sub- 
sonic models show no correlation between column den- 
sity and magnetic field. Supersonic cases are able to 
reach higher field strength for higher densities due to 
the entanglement of the field lines with density clumps, 
which is similar to both the 3D correlation of magentic 
energy vs. p in Figure [3] and the 2D correlation in Fig- 
ure [SI The difference between Figures [5] and [6] is the 

^ The magnetic field _L to the line of sight can be estimated with 
the Chandrasechar-Fermi technique. 
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Fig. 4. — The 3D correlation of the local Alfvenic Mach number (Ma) vs. normalized density for a log log scale. The first row consists of 
globally super-Alfvenic cases while the bottom row is globally sub-Alfvenic. The mean Alfvenic number is always nearly the global value 
for the model. Images are ordered left to right as supersonic to sub-sonic. Blue contours indicate regions of high data counts while red and 
yellow have lower counts. We show the actual values of Ma for ease of reading. 



orientation of the field with respect to the plane of pro- 
jected column density. When the field is perpendicular 
to this plane we see a higher degree of linearity then in 
the parallel case. This is because density clumps have 
more freedom when the compressions are perpendicular 
to the field lines and thus clumps are able to better trap 
magentic energy. When the field is parallel to the plane, 
plasmas tend to fiow along field lines and cannot trap 
them. Correlations are stronger for even subsonic cases 
when the column density is perpendicular to the field and 
the correlation coefficients of Table [6] confirm this. We 
must keep in mind that while the field is oriented par- 
allel and perpendicular to the column density integrated 
along the x direction for Figures [5] and [6] respectively, the 
mean volume magnetic field is parallel to the x direction 
in our simulations. 

Figure [71 shows the correlation of column density along 
z-direction and the magnetic field component parallel to 
the line of sight. Figure [8] shows the correlation of column 
density along z-direction and the magnetic field compo- 
nent perpendicular to the line of sight. Both correlations 
are obtained along the z direction, i.e. perpendicular 
to Bq. Figure [71 is similar to Figure [5] for models with 
Ma = 2.0. and Figure [8] is similar to and shows steeper 
correlations then [6] for all models. This occurs because 
the line of sight is chosen to be ± Bq and the observed 
magnetic field is simply the random/perturbed compo- 
nent. We must keep in mind that the density struc- 
tures in MHD turbulence are, in general, filamentary and 
anisotropic regarding the magnetic field orientation. 

As a consequence, column density maps along x and z- 
directions will be different. This effect is more prominent 
on sub-Alfvenic models, as the magnetic field lines are 
slightly changed due to turbulent motions. We visually 



see that a strong magnetic field compresses matter more 
when densities are perpendicular. This is clearly shown 
since the plots in Figure [8] are the most linear of all 
four B vs. column density correlations. Also, despite the 
orientation of the field, it is noticeable from Figures [5][8] 
that supersonic models present smaller dispersion over 
the column density-magnetic field space. 

In order to compare the observational measures of ve- 
locities from spectral lines with the synthetic velocity dis- 
persions, we studied the correlations of column density 
and velocity dispersion parallel and perpendicular to the 
line of sight. We present the figures for dispersion per- 
pendicular to the LOS (shown in Figure [9]). We calculate 
velocity dispersion by integrating ayp along the LOS. Su- 
personic models present stronger correlation along the z 
LOS then along the x LOS due densities being perpen- 
dicular to the mean magentic field. The effect of this 
orientation is explained in Figures [3 and [H For cases 
parallel to the magnetic field, matter is confined by the 
magnetic field and the fiows tend to be along it, thus 
there is ultimately less dispersion seen in these cases. 
When column density is perpendicular to the field, as 
shown in Figure [9l dispersion increases with density for 
supersonic cases. The higher velocities give rise to higher 
compression of matter due to shock waves. This is evi- 
dent upon comparison between supersonic and subsonic 
models Exactly how to relate the statistics of observed 
velocity line profiles to simulations is not straight for- 
ward (see Esquivel & Lazarian 2005). The interaction 
between density fiuctuations and velocity in synthetic 
MHD cubes must first be characterized before more com- 
plicated statistics can be applied and comparisons with 
observations are made. We show the correlation between 
observable velocity dispersions in Figure [TOl We study 
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Fig. 5. — The 2D correlation of the integrated magnetic field component parallel to the line of sight vs. column density in the x direction. 
The first row consists of super- Alfvenic models while the bottom row is sub-Alfvenic. Images are ordered left to right as supersonic to 
subsonic. Blue contours indicate regions of high data counts while red and yellow have lower data counts. 



how the actual dispersion of velocity correlates with the 
measure available through spectral observations of the 
Doppler shifted lines. We understand the dispersion of 
velocity along the z axis to be Vz dispersion= (Jy^p^ - Due 
to shock waves, we expect strong density fluctuations to 
be present for supersonic cases. This will result in gener- 
ally nonlinear correlations of velocity dispersions for con- 
stant vs. non constant p. However, for subsonic cases, 
density fluctuations are less prevalent, and linear correla- 
tions should be present. Considering dispersion, models 
with higher Mach numbers should reach higher values 
of velocity dispersion, by deflnition. All of these trends 
are clearly shown in Figure [TOl Supersonic models show 
no correlation due to density fluctuations while subsonic 
models are linear. The stronger the magnetic fleld, the 
tighter the correlation. This trend is conflrmed in Figure 

El 

4.3. Correlation Coefficients 

Finally, we present a table of select log-log correlation 
coefficients with error bars in Tables [5] and [6] to see 
quantitatively the degree of correlation in our plots. Ul- 
timately, we seek to set an order on which models show 
the highest correlations. These correlation coefficients 
tell us the degree of the polynomial equation fitted to 
quantities plotted. In order to obtain these coefficients 
we computed the linear least-squared fit of the correla- 
tions, with error bars determined from the 1-sigma un- 
certainty estimates. For magnetic energy vs. density, 
highly supersonic cases show positive correlations while 
subsonic present anti-correlation. However, it is of in- 
terest to note that the cases with the largest correlation 
differ in sonic Mach number, but both have A^^=2.0. 
We see that from our results that both the Ms and Ma 
play a roll in the correlation of B and p, due to shocks cre- 
ating density clumps in supersonic cases. The A1a=2.0, 



A^A=0.7 case shows the closest to Alfvenic perturba- 
tions {B ^ \^p)' In general, magnetic energy cannot be 
correlated with p with a simple polynomial expression. 
For kinetic energy vs. density, all show negative correla- 
tions with supersonic cases being correlated the greatest, 
showing that no strong correlation trend was obtained. 
For magnetic field parallel to LOS vs. x column density 
(Sa^), the table shows that the most important factors in 
obtaining a strong correlation is the sonic Mach number 
and the strength of the field. Supersonic cases with large 
magnetic field are generally linear while supersonic cases 
with small magnetic field show an trend. Looking at 
magnetic field perpendicular to LOS, it is apparent all 
models are near linear. The orientation of the field per- 
pendicular to the LOS creates correlations that are more 
nearly linear and this is evident in both the plots and the 
coefficients of B± vs. S^. We also show the correlation 
coefficients for Ma vs. p. No polynomial relationship is 
seen for any model. For velocity dispersion and S^, a 
strong linear relationship is established which is not seen 
for the models that are along T^x (parallel to the mag- 
netic field). It is evident that the orientation of the field 
with respect to the LOS is very important in obtaining 
a strong correlation with dispersion. 

5. BISPECTRUM 

5.1. Definition and Calculation of bispectrum 

As turbulent vortexes evolve, they transfer energy from 
large to small scales. In this case, wave- wave interactions 
generate a hierarchical turbulent cascade as ki-\-k2 ks. 
For incompressible fiows, under Kolmogorov's assump- 
tions, we have ki k2 = k and ks — 2k. For com- 
pressible and magnetized fiows, this becomes more com- 
plicated and nonlinear wave-wave interactions may take 
place, mainly for ki ^ k2- Also, we expect MHD turbu- 
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Fig. 6. — The 2D correlation of the integrated magnetic field component perpendicular to the line of sight vs. column density in the 
X direction. The first row consists of super-Alfvenic models while the bottom row is sub-Alfvenic. Images are ordered left to right as 
supersonic to subsonic. Blue contours indicate regions of high data count while red and yellow have lower counts. 
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Fig. 7. — The 2D correlation of the integrated magnetic field component parallel to the line of sight vs. column density in the z direction. 
The first row consists of super-Alfvenic models while the bottom row is sub-Alfvenic. Images are ordered left to right as supersonic to 
subsonic. Blue contours indicate regions of high data counts while red and yellow have lower counts. 
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Fig. 8. — The 2D correlation of the integrated magnetic field component perpendicular to the line of sight vs. column density in the 
z direction. The first row consists of super-Alfvenic models while the bottom row is sub-Alfvenic. Images are ordered left to right as 
supersonic to subsonic. Blue contours indicate regions of high data counts while red and yellow have lower counts. 
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Fig. 9. — The 2D correlation of velocity dispersion parallel to the line of sight (i.e. along the z direction) vs. column density in the z 
direction. The mean magentic field is 1. to LOS. The first row consists of super-Alfvenic models while the bottom row is sub-Alfvenic. 
Images are ordered left to right as supersonic to subsonic. Blue contours indicate regions of high data counts while red and yellow have 
lower counts. 
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Fig. 10. — The 2D correlation of dispersion velocity for constant p vs. the dispersion velocity along the Z direction (i.e. mean magnetic 
field -L to the LOS). The first row consists of super- Alfvenic models while the bottom row is sub-Alfvenic. Images are ordered left to right 
as supersonic to subsonic. Blue contours indicate regions of high data counts while red and yellow have lower counts. Supersonic cases 
show evidence of nonlinear density fiuct nations while subsonic cases are linear. 



TABLE 5 

Table of log-log Correlation Coefficients 



Model 


Magnetic Energy vs. p 


Kinetic Energy vs. p 


B\\ vs. Ex 


B± vs. Ex 


1 


-1.1700 ± 1x10-4 


-0.5810 ±3x10-4 


0.020 ± 0.003 


1.200 ±0.004 


2 


-0.10000 ^4x10-5 


-0.1500 ± 1x10-4 


0.950 ±0.001 


1.240 ±0.002 


3 


0.00200 ± 2x10-5 


-0.05000 ±8x10-5 


0.990 ±0.001 


1.240 ± 0.002 


4 


-1.0600 ±4x10-4 


-0.6330 ± 3x10-4 


0.720 ± 0.009 


0.920 ±0.003 


5 


0.1500 ±1x10-4 


-0.1300 ±1x10-4 


2.270 ± 0.004 


1.310 ±0.001 


6 


0.3400 ± 8x10-5 


-0.02000 ± 7x10-5 


2.110 ±0.003 


1.350 ±0.001 



TABLE 6 

Table of log-log Correlation Coefficients 



Model 


B\\ vs. E^ 


B^ vs. Ez 


Ma vs.p 


^dispersion VS. E^ 


^dispersion Ej\" 


1 


-0.13 ±0.01 


1.3600 ± 2x10-4 


0.7980 ± 2x10-4 


0.810 ±0.006 


0.350 ±0.007 


2 


1.550 ±0.004 


0.9300 ± 6x10-4 


0.47300 ± 6x10-5 


1.180 ±0.003 


-0.170 ±0.004 


3 


1.570 ±0.003 


0.9800 ± 5x10-4 


0.47100 ±4x10-5 


1.190 ±0.002 


-0.120 ±0.004 


4 


2.130 ±0.008 


1.300 ±0.003 


0.7150 ±3x10-4 


1.19 ±0.005 


0.620 ± 0.005 


5 


2.330 ± 0.004 


1.320 ±0.001 


0.35700 ±9x10-5 


1.1200 ±2.7x10-4 


0.230 ± 0.003 


6 


2.020 ± 0.003 


1.360 ±0.001 


0.32020 ±6x10-5 


1.270 ±0.002 


0.616 ±0.005 



lence to present more wave modes than purely hydrody- 
namical flows and, in this case, the energy cascade can 
be much more complicated. In order to study this phe- 
nomenon and characterize both amplitude and phase in 
the turbulent cascade, the analysis of the three-point cor- 
relation function, or bispectr um, is required (see Barnett 
I2QQ2I : iMasahiro fc Bhuvneshi r2QQ4). 

The bispectrum technique characterizes and searches 
for nonlinear interactions, which makes it a timely study 
for interstellar MHD turbulence. However, unlike the 



standard power spectrum, the bispectrum is a complex 
function and carries phase information. In order to com- 
pute either the power spectrum or bispectrum, one must 
first compute the Fourier transform of the signal. In the 
case of either a discrete or continuous Fourier transform 
we are left with information that includes both the mag- 
nitude and the phase of the signal at a given frequency. 
The power spectrum measures solely the magnitude as 
a function of the frequency and discards the phase in- 
formation. The power spectrum gives information about 
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power distributions over all frequencies and looses infor- 
mation about wave- wave interactions and phase. Signals 
containing multiple frequencies and phases are better de- 
scribed by the bispectrum. 

The bispectrum is closely related to the power spec- 
trum. The Fourier transform of the second-order cu- 
mulant, i.e. the autocorrelation function, is the power 
spectrum while the Fourier transform of the third order 
cumulant is known as the bispectrum. In a discrete sys- 
tem, the power spectrum is defined as: 

P{k)= Mk)-A*{k) (6) 

k=const. 

In a similar way, the bispectrum can be defined as: 

k\=const k2=const 

(7) 

where ki and k2 are the wave numbers of two interacting 
waves, and A{k) is the original discrete time series data 
with finite number of elements with A^{k) representing 

the complex conjugate of A{k). As is shown in Equa- 
tion [71 , the bispectrum is a complex quantity which will 
measure both phase and magnitude information between 
different wave modes. 

In practice, our calculation of the bispectrum involves 
performing a Fast Fourier transform of density or column 
density and the application of Equation [3 Although we 
are primarily interested here in density fiuct nations, the 
bispectrum can also be calculated for other fields as well. 
We randomly choose wavevectors and their directions, ki 
and /c2 and iterate over them, calculating /cs, which de- 
pends on ki and k2 since ki -\- k2 = ks. We limit the 
maximum length of the wave vectors to half of the box 
size. We normalize direction vectors to unity, calculate 
positions in Fourier space, and finally, compute the bis- 
pectrum, which yields a complex number. We then aver- 
age this bispectrum over all frequencies and plot bispec- 
tral amplitudes. This gives us information on the degree 
of mode correlation in the system. 

In this section we will explore the bispectrum of den- 
sity and column density for the data cubes. Because this 
is the first application of the bispectrum to ISM related 
studies and we would like to have a comparison with 
observations, exploring fields other then density and col- 
umn density is beyond the scope of this paper. Although 
the application of the bispectrum to velocity and mag- 
netic fields would be interesting and provide a beneficial 
follow up to densities, we focus here only on density and 
column density maps in order to make our study com- 
parable with observational column densities. Bispectral 
analysis can be applied to observational data in a similar 
way as the synthetic data presented in this paper utiliz- 
ing Equation [3 Our model does not include self-gravity 
or external gravity, and only considers the interplay be- 
tween the gas pressure and magnetic pressure of isother- 
mal gas. By examining the bispectrum of density and 
comparing it with column densities, we can characterize 
turbulent fiows on resolutions that are realistically seen 
in observational data. We will look at different mod- 
els of turbulence including a Gaussian, pure hydro, and 
different MHD regimes. We will examine differences be- 



tween the bispectrum of density and column density and 
discuss what information about the sonic Mach number 
and magnetic field that might be gained by applying the 
bispectrum to observable data. 

5.2. Bispectra of Density and Column Density 

We ran the bispectrum analysis for the last snapshots 
of the density and column density cubes on x and y- 
directions for all MHD models. We also performed the 
bispectrum calculation for two comparison models: a su- 
personic hydrodynamical simulation and a random syn- 
thetic Gaussian distribution of density. The density and 
column density integrated along x and y bispectra are 
shown in Figure [Til the left, center and right columns, 
respectively. In the the first row we present the results 
obtained for the Gaussian distribution, followed by the 
hydrodynamical case and then the MHD models. The 
models are labeled for the different turbulent regimes 
with varying values of and Mg- We plotted con- 
tours and their corresponding color amplitudes in order 
to see both differences in the shape and magnitude of the 
bispectrum. 

Upon examination of Figure [TT] one sees a general cor- 
respondence between the bispectrum of the column den- 
sity and the underlying 3D bispectrum. All cases show a 
prominent amplitude of the bispectrum and thus a high 
correlation of modes, at ki = k2. This is to be expected 
since this is a trivial case where wave numbers will always 
show correlation. However, the amplitudes at ki ^ k2 are 
different for each turbulent regime. The more circularly 
shaped the isocontours and stronger the amplitudes are, 
the more highly correlated the modes are. The Gaus- 
sian distribution shows almost no correlation (very weak 
amplitudes) for k\ ^ k2 and the distribution is mostly 
localized in the diagonal line (i.e. where ki = /C2). This is 
because the distribution of density is entirely random and 
thus modes should lack association. The supersonic hy- 
dro model is similar to a super- Alfvenic supersonic MHD 
model, as expected, and is more strongly correlated then 
the Gaussian density distribution. The MHD models also 
present broader distributions over the /ci, k2 plane, thus 
showing high correlation. 

The bispectrum of MHD density is shown in the first 
column starting at the third map down. Larger wave 
numbers (representing larger frequencies and smaller 
wavelengths) are less correlated in all models, since 
higher wave numbers are present near the end of the 
energy cascade when more of the systems energy has dis- 
sipated. For the models with the same Ma but varying 
Ms^ we found an increase in the wave- wave coupling 
under the supersonic regime. This phenomenon is ex- 
pected, as interactions between different modes generate 
small scale shocks. Shocks in supersonic models greatly 
increase wave- wave correlations over subsonic models, 
which are incompressible. The supersonic models of den- 
sity (1st column, last and second to last rows) show such 
strong correlations that the strength of the ki 7^ /c2 is 
comparable to the ki = k2 line. Interestingly, comparing 
models with same sonic Mach number, the sub- Alfvenic 
cases show increased wave- wave coupling for ki 7^ /c2. 
Here, non-linear coupling is increased by the magnetic 
field. Comparing density models with Ms = 7.0 we see 
that a stronger magnetic field shows slightly higher am- 
plitudes and has a more circular isocontour profile. Com- 
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Fig. 11. — Above is the contour analysis of the bispectrum for density and column density. These figures show the degree of 
correlation between wavenumbers ki and k2. The first column shows density, the second shows column density parallel to the 
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paring Ms = 0.7 models, it is easy to see that stronger 
amplitudes exist for the cases with Ma = 0.7. Looking 
specifically at the second and third columns for the su- 
personic models, the amplitudes of the Ma = 0.7 models 
are much stronger then that of the Ma = 2.0 models. 
It is clear from these models that a stronger magnetic 
field causes enhanced correlation of modes. For super- 
Alfvenic models, the large scale (low k) magnetic field 
configuration is destroyed and the MHD modes operate 
at smaller scales. Both effects reveal that the energy 
cascade, generated or amplified by ki ^ k2 interactions, 
may work differently on the turbulent regimes, and Kol- 
mogorov scalings may not well-characterize them. 

Regarding observable parameters, the column density 
bispectra do not show such prominent differences com- 
pared to the density analysis. However, the column den- 
sity bispectra present larger noise because of worse statis- 
tics (512^ maps compared to the 512^ density cube). 
However, the supersonic and sub-Alfvenic model shows 
a notably different distribution, with larger amplitudes 
for /ci 7^ , similarly to the density bispectrum. Smaller 
amplification is seen in the subsonic cases, particularly 
the subsonic super- Alfvenic model. It is clear from the 
bispectra of density and column density that supersonic 
models present the highest correlation of modes while 
models with a high magnetic field also enhances mode 
correlation. If this is true for synthetic data then the 
bispectrum technique could be used to characterize mag- 
netized turbulence in the ISM. Interestingly, the broad- 
ening of bispectrum distribution by the magnetic field is 
mostly independent on the orientation of the magnetic 
field lines regarding the line of sight, however differences 
are clear for supersonic sub-Alfvenic cases. The bispec- 
trum is slightly broader for the x direction since the mean 
magnetic field is along this axis. The column density 
bispectra present larger noise because of worse statistics 
(512^ maps compared to the 512^ density cube). From 
our studies it is obvious that the bispectrum requires very 
high resolution data in order to yield meaningful results. 
Noise and other outside infiuences can completely mask 
a desired signal information in the bispectrum. Noise 
will produce a 'false' bispectral map, showing higher cor- 
relations for larger wave numbers. Therefore, bispectral 
analysis of observational maps can be used to reveal the 
turbulent regimes operating in molecular clouds, but only 
if high resolution maps with minimal noise are available. 

5.3. Bispectra and spectra 

Quantitative studies of compressible MHD turbulence 
in Cho & Lazarian (2002, 2003), concentrated on decom- 
position of turbulence perturbations into slow, fast and 
Alfvenic modes and studies of their spectra. Such a de- 
composition requires the knowledge of the local magnetic 
field, which is not the case for the data averaged along 
the line of sight. Therefore from observations the to- 
tal spectra of turbulent perturbations are available. The 
density spectrum is readily available from column density 
measurements, while the velocity spectrum requires more 
sophisticated techniques, like Velocity Channel Analysis 
(VGA) or Velocity Coordinate Spectrum (VCS) to be 
used (see Lazarian & Pogosyan 2000, 2006, 2008). Dif- 
ferent numerical studies show that for MHD turbulence, 
the velocity spectrum gets steeper for velocity and shal- 



low for density p as the Ms increases^^. KLB showed 
that the dependences of the density spectrum on A^^ is 
less prominent than on Ms- 

Above we studied the bispectrum of density and no- 
ticed the dependence of the measure of Ms- We see that 
the variations of the bispectrum with Ma for all ki and 
/c2, and find that the phases of the turbulent fiuctuations 
are more correlated in the presence of the dynamically 
important magnetic fields. 

6. DISCUSSION 

6.1. Comparison to KLB and other studies. 

KLB presented an extensive statistical analysis of 
MHD density for cubes of resolution 128^, 256^, and 
512^. However, they only examined the skewness and 
kurtosis for the lower resolution model (128^). In this 
work we study cubes of 512^ resolution in order to com- 
pare with the statistics of lower resolution models. KLB 
examined a larger range of Ms values and noticed that 
the asymmetry grows with Ms, but only if Ms > 0.5. 
In this work, we study models with 0.7 < Ms < 7.0 
and we find that statistical moments such as skewness 
and kurtosis are easy to obtain and are very useful for 
characterizing sonic Mach number in a magnetized tur- 
bulent system. We confirm the results of KLB in find- 
ing that the moments all have strong dependence on the 
sonic Mach number. As turbulence becomes increasingly 
supersonic, the skewness of density and column density 
increases, which is agreement with KLB. We also see a 
very rapid growth in kurtosis of density and column den- 
sity with sonic Mach number. Both column density and 
density show increasing asymmetric distributions with 
increasing Ms- If one observes a given skew and kurto- 
sis of column density, the value of A^^ can be inferred 
and related to density, which could be very useful when 
applied to observational data. It seems that in both our 
study and in KLB, Ma does not affect the moments of 
density. 

Our study of magnetic fields and column density corre- 
lations in isothermal nonlinear MHD waves has similar- 
ities to that of Passot & Vazquez-Semadeni(2003) who 
found that the relation of magnetic field strength B on 
density p depends on whether the MHD wave is a slow or 
fast mode. Their results confirmed that, except for cases 
of insufficient field fiuctuations, the magnetic pressure 
behaves as: 5^ c± ci — for slow modes (low Ma) and 
B'^ ^ p^ for fast modes. However, while this study in- 
corporated a large variety of values for the Alfven Mach 
number (0.073 < Ma < 7.29), they lacked sampling 
of sonic Mach numbers, only focusing their analysis on 
supersonic values, 2.72 < Ms < 7.28. They also only 
investigated 2D models, while in this study we use 3D 
simulations of magnetic energy and density as well as 
2D correlations of magnetic field and column density at 
higher resolutions. In this study we also explore sub- 
sonic as well as supersonic cases and find that the sonic 
Mach number is a critical parameter in determining how 
B and p correlate, since it determines whether or not 
density clumps will develop. In agreement with Passot 
& Vazquez-Semadeni(2003), we also find that a larger 

As a result the quantity p^l^v stays invariant to M.s both in 
hydro (Kritsuk et al. 2006) and MHD (Kowal & Lazarian 2006) 
turbulence. 
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value of Ma yields a tighter correlation between density 
and magnetic pressure, especially at higher density. This 
effect can easily be seen by comparing the 3D correlations 
in Figure [3l the Ads = 7.0 cases. However, Figures [Ml] 
confirm a more linear relationship for slow modes ( i.e. 
low Ma) over larger Ma^ which contain both slow and 
fast modes. We also see that the orientation of the field 
along the LOS can effect these results and we confirm 
Passot & Vazquez-Semadeni's result that the angle be- 
tween the magnetic field and the wave propagation will 
play a role in determining the importance of the modes. 
Subsonic models generally show little correlations, yet 
in the 2D correlations we show that when the LOS is 
perpendicular to the magnetic field and mean field (Fig- 
ure [8|), subsonic correlations begin to show a degree of 
linearity. 

While higher order moments show virtually no depen- 
dence on Ma^ correlations, on the other hand, can yield 
information regarding the dynamical importance of the 
magnetic field in the evolution of densities. According 
to Padoan and Norlund(1999), cloud dynamics only be- 
come strongly affected by magnetic field in very dense 
regions and on small scales. It was also further stated 
in Padoan et al.(2004) that, while the volume averaged 
magnetic field strength in a molecular cloud has never di- 
rectly been measured, there is evidence from the power 
spectrum of gas densities that super- Alfvenic turbulence 
is dominate in these star forming regions. From the 3D 
correlation in Figure [3] it is clear that magnetic fields play 
an important role in shaping density. This is especially 
true of supersonic models where a correlation of mag- 
netic energy and density exists. Higher density regions 
are good candidates for star formation, thus it may seem 
promising that higher magnetic field strength could in- 
dicate more star forming regions. Also, from Figure [H 
it is clear that even though turbulence may be globally 
super or sub- Alfvenic, dense dumps influence magnetic 
field to a point where the local Alfven Mach number can 
change greatly. In the case of super- Alfvenic turbulence, 
high density results in smaller values Ma (although they 
generally remain super- Alfvenic) while in sub-Alfvenic 
turbulence. Ma steadily increases with density to the 
point where turbulence becomes super-Alfvenic. How- 
ever, even though clumps in globally sub-Alfvenic clouds 
are able to become locally super-Alfvenic, this does not 
mean that the cloud as a whole becomes super-Alfvenic. 
We do not see any evidence that the cloud as a whole is 
unaffected by a strong magnetic field, although we do see 
that indeed a magentic field has a strong impact on small 
scale turbulence. Therefore, our globally super-Alfvenic 
models support Padoan and Norlund(1999)'s claim that 
magnetic field plays a leading role in the dynamics of 
very dense clumpy regions (i.e. as density increases, 
MAdecreases) , however there is no evidence for this in 
our sub-Alfvenic models, which show that on small scales 
density clumps are increasingly super-Alfvenic, and thus 
have a lower B field. Overall, these correlations only give 
us information on how the magnetic field effects local 
density clumps and not whole cloud. However, it is per- 
haps promising that our results tend to be locally super- 
Alfvenic for high densities, despite the global magnetic 
field. More studies will be needed in order to confirm 
Padoan and Norlund(1999) results. 

While density statistics are improtant to understand- 



ing turbulent processes, the observational velocity infor- 
mation provided by radio spectroscopic data should be 
utilized in order to gain a broader study of turbulent in- 
teractions. The difficulty in directly comparing synthetic 
velocity dispersions with observational spectral velocities 
lies in the fact that the emissivity of a spectral line has 
dependence on velocity and density. In order to bet- 
ter understand this, we calculate correlations of velocity 
dispersion vs. column density and also investigate how 
dispersion of velocity correlates with the measures avail- 
able through Doppler shifted spectral lines. Lazarian 
et al.(2001) studied emissivity statistics and the density- 
velocity correlation to better understand how density and 
velocity interact in a turbulent media. In this study, 
we confirm that a moderate correlation between column 
density and velocity exists, specifically seeing that higher 
velocities indicate higher compression of matter. How- 
ever, we further find here that these correlations gain 
strength when the LOS is ± to the magnetic field. The 
Lazarian et al. study used lower resolution cubes and 
only examined 3D cases, while here we investigate 2D 
cases at higher resolutions. In general, further studies 
must be made in order to understand how best to relate 
observational velocity line profiles to simulated velocity 
dispersions. 

6.2. Density Correlations and Relation to Observables: 
2D vs. 3D 

In order to compare correlations of 3D densities to 
quantities more similar to what is available through ob- 
servations (i.e. column density, velocity dispersions), we 
provide 2D column density correlations. We do indeed 
see several similarities between 3D correlation maps and 
their 2D counterparts. For instance, when comparing the 
correlation of density and magnetic energy (Figure [3]) 
with correlations of magnetic field and column density 
(Figures [SJIH]), we see that for supersonic cases, magnetic 
energy {B'^) and magnetic field both increase with den- 
sity and column density. This implies that for 2D and 
3D cases, magnetic fields become trapped in high den- 
sity clumps. However, because Figure [3] is a volume 
plot and we plot magnetic energy instead of field, we see 
differences. Most notably, we see a 5 ^ p for 2D maps 
and increases more exponentially with p for super- 
sonic cases. We can see this trend from Tabled with 2D 
correlations being linear with cases with the field per- 
pendicular to the LOS being the most linear. Therefore, 
although the magnetic field and energy both see a similar 
growth with density for supersonic turbulence, the rate 
of growth is different. Subsonic cases for both 3D and 
2D correlations both showed the least correlation. 

Specific kinetic energy {£k) vs. density showed no spe- 
cial relationship. However, the kinetic energy weighted 
by density (Ek) shows a nearly linear correlation for su- 
personic cases. We see a similar behavior for velocity dis- 
persions vs. column density. These dispersions are also 
weighted by density, since observers must also account 
for the column density when a measure of velocity disper- 
sion is made. Because of density clumps, both 2D and 3D 
quantities increase with dispersion and kinetic energy, re- 
spectively. For both correlations, subsonic models show 
less linear correlations. However, despite many similari- 
ties, it is not necessarily straightforward to compare the 
two, since the velocity dispersion seems to be effected by 
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the orientation of the magnetic field with the LOS. We 
have not tested this effect on kinetic energy density. It is 
obvious from the 2D correlations that the field orienta- 
tion effects how correlated dispersions are with column 
density. When the field is perpendicular to the LOS, the 
most linear correlations are obtained, as is evident from 
TableO Dense clumps can drag the field lines and, there- 
fore, present a more isotropic velocity field. This effect 
results in an increase in the dispersion of the velocity 
measured along the LOS. 

6.3. 2D and 3D Bispectrum 

The bispectrum can give us valuable informa- 
tion regarding how the modes of a nonlinear sys- 
tem correlate. The power spectrum has been 
used in nun ierous papers to analyze o b serva- 
tions, fsee lArmstrong. Rickett. fc Spanglerl Il995l: 
Stanimirovic et aTl 1 19991 : IStanimirovic fc LazarianI 
2QQ1D . In Section E we analyzed the bispectrum for 
densities and column densities and we have found several 
interesting features which may provide new ways of 
analyzing data. Different turbulent parameters display 
very different bispectra, and thus the bispectrum maybe 
be used to characterize Mach number and magnetic field 
in observational data. 

For instance, the bispectrum of density and column 
density gives information as to how shocks and mag- 
netic fields affect turbulence. It has been shown by KLB 
and Beresnyak et al.(2005) that in supersonic turbulence, 
shocks produce compressed density and a shallower spec- 
trum. Looking at Figure [TT] for density (first column) it is 
clear that these shocks play a crucial role in the correla- 
tion of modes. The subsonic cases show little correlation 
between any points except the case of ki = /c2 . The com- 
pressed densities from supersonic turbulence are critical 
for correlations between wave numbers, since waves will 
be compressed together and therefore have more inter- 
actions. Also interesting to consider is the differences in 
sub-Alfvenic and super-Alfvenic cases for density. We 
find excellent agreement with KLB in that density struc- 
tures are well correlated to the presence of magnetic 
fields. For super-Alfvenic cases, it is clear from Fig- 
ure [TT] that the system lacks the stronger correlations 
that are a characteristic of the sub-Alfvenic models. For 
super-Alfvenic simulations, correlations in frequency are 
not as readily made due to large dispersion of density 
structure. The bispectrum of supersonic hydro models 
is similar to the super-Alfvenic, supersonic cases. It is 
clear from these cases that the presence of a magnetic 
field assists in the correlation of modes. 

In order to relate to observations, the bispectrum of 
column densities is provided in order to compare with 
density. It should be noted that while our model does 
not include self or external g:ravity, it has been shown 
by (Elmegree n fc Sca'lQll2QQ4l ) that self-gravity partitions 
the gas into clouds which contributes to scale-free mo- 
tions generated by turbulence. Thus, self-gravity en- 
hances small scale structure and our results should be 
relatively unaffected, as we are primarily concerned with 
large scale turbulent cascades. For both density and col- 
umn density, strongest correlations are found in the case 
of supersonic sub-Alfvenic turbulence. Differences in im- 
age quality between density and column density maps 
arise from the fact that density cubes have higher reso- 



lution {N^ vs. A/"^). Because the magnetic field enhances 
correlations, the bispectrum could be used to character- 
ize the magnetic field in studies similar to Goodman et 
al.(1995); Padoan fc Nordlund (1999). 

The bispectrum has proved to be an important addi- 
tion to the tools of statistical studies of MHD turbu- 
lence. This new tool could be used to characterize the 
Mach number of gas in the interstellar medium as well 
as the magnetic field. However, it should be noted that 
very high quality data is required since the multipoint 
statistics are known to increase noise. For the interstel- 
lar medium it has been suggested in ([LazarianI 1 1999) to 
compare various regions of the sky using the bispectrum 
to search for like signals. These studies are to be done 
elsewhere. 

7. CONCLUSIONS 

In this paper we investigated several different statistics 
of density structure of compressible MHD turbulence. 
We examined the skewness and kurtosis of the data as 
well as studied several different measures of correlations. 
In particular, we examined the bispectrum of compress- 
ible MHD turbulence. We found: 

• We confirm KLB's results with higher resolution 
cubes that sonic Mach number has dependence on 
the 3rd and 4th order moments. 

• Correlations of kinetic energy density vs. density 
are linear for supersonic cases and show no corre- 
lation for subsonic cases 

• Magnetic energy increases with density clumps for 
supersonic turbulence. 

• For the super-Alfvenic cases, the local Alfven Mach 
number, i.e. Ma = of clumps increases up to 
the mean density then falls of as B increases with 
increasing density clumps. 

• For sub-Alfvenic cases, the local Ma increases with 
density for both supersonic and subsonic cases. 

• Correlations of magnetic field vs. column den- 
sity are linear for supersonic turbulence and are 
strongest when the field is perpendicular to the 
LOS. 

• Correlations of velocity dispersion vs. column den- 
sity show a degree of linearity and are strongest 
when the mean magnetic field is perpendicular to 
the LOS. 

• Subsonic velocity dispersions are not strongly ef- 
fected by density fiuctuations. 

• Applying the bispectrum to density fields of MHD 
turbulence we find that: 

1. There are strong correlations for cases where 
ki = k2 

2. There are virtually no correlations for ki ^ k^ 
for subsonic cases. 

3. Supersonic models show the strongest k\ ^ k2 
correlations. 
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4. The introduction of a magnetic field enhances 
correlations for all ki and . 



B.B. is supported by the NSF funded Research Expe- 
rience for Undergraduates (REU) program through NSF 



award AST-0453442. B.B. is thankful for valuable dis- 
cussions with Dr. James T. Lauroesch and Thiem Hoang. 
D.F.G. thanks the Brazilian agencies FAFESF (No. 
06/57824-1) and CAFES (No. 4141067). We also thank 
the financial support of the NSF (No. AST0307869) and 
the Center for Magnetic Self- Organization in Astrophys- 
ical and Laboratory Flasmas. 



REFERENCES 



Armstrong, J.W., Rickett, B.J. & Spangler, S.R., 1995, ApJ, 

443,209 

Ballesteros-Paredes, J., Klessen, R. S., Mac Low, M.-M., &; 

Vazquez- Semadeni, E. 2007, in Protostars and Planets V, ed. 

B. Reipurth, D. Jewitt, &; K. Keil, 6380 
Beresnyak, A., Lazarian, A. & Cho, J., 2005, ApJ, 624, L93 
Barnett, A. G. 2002, PhD Thesis, School of Physical Sciences, The 

University of Queensland. 
Cho, J. &; Lazarian, A. 2002, Phys. Rev. Lett., 88, 5001 
Cho, J. &; Lazarian, A. 2003, MNRAS, 345, 325 
Cho, J., Vishniac, E., Beresnyak, A., Lazarian, A., Rye, D., ApJ 

submitted 

Draine, B.T., & Lazarian, A. 1998a, ApJL, 494, L19 
Elmegreen, B. & Scalo, J. 2004, ARA&A, 42, 211 
Esquivel, A., Lazarian, A., Pogosyan, D., & Cho, J. 2003, MNRAS, 
342, 325 

Esquivel A., Lazarian,A., 2005, ApJ, 631, 320 
Evans, N. J. 1999, ARA&A, 37, 311 

Falceta-Gongalves, D., Lazarian, A., Kowal, G. 2008, ApJ, 679, 537 
Falceta-Gongalves, D., de Juli, M. C. & Jatenco-Pereira, V. 2003, 
ApJ, 597, 970 

Fry, J. N. 1998, Annals of the New York Academy of Sciences 848, 
62 

Fry, J. N. & Seldern, M. 1982, ApJ, 259, 474 

Gerakines, P. A., Whittet, D. C. B. & Lazarian, A. 1995, ApJ, 455, 
171 

Girart, J. M., Rao, R. & Marrone, D. P. 2006, Science, 313, 812 
Goodman, A. A., Jones, T. J., Lada, E. A. & Myers, P. C. 1995, 

ApJ, 448, 748 
Haffner, L. M. et al. 2003, ApJS, 149, 405 

Hih, A. S., Benjamin, R. A., Reynolds, R. J. & Haffner, L. M. 2005, 
BAAS, 1302 

Hill, A. S., Benjamin, R. A., Kowal, G., Reynolds, R. J., Haffner,. 

L. M. & Lazarian, A., 2008, ApJS 
Hiltner, W. A. 1956, ApJS, 2, 389 

Intrator, T., Meassick, S., Browning, J., Majeski, R. &; 
Hershkowitz, N. 1989, Phys. Fluids B, 1, 271-273 

Kowal, G., Lazarian, A. & Beresnyak, A., 2007, ApJ, 658, 423 

Lazarian, A. 1999, Plasma Turbulence and Energetic Particles 
in Astrophysics, Proceedings of the International Conference, 
Krakow (Poland), 5-10 September 1999, Eds.: Michal Ostrowski, 
Reinhard Schlickeiser, Obserwatorium Astronomiczne, 
Uniwersytet Jagielloznski, Krakow 1999, p. 28-47., 28 

Lazarian, A. & Cho, J. 2004, "Basic Properties of Compressible 
MHD: Implications for Molecular Clouds", Astrophysics and 
Space Science, 292, 29-43 

Lazarian, A., Kowal, G., Beresnyak, A., 2008, ASP Conf. Ser. Vol 
385 p.3 

Lazarian, A., Pogosyan, D., & Esquivel, A. 2002, ASP Conf. Proc. 
276, Seeing Through the Dust: The Detection of HI and the 
Exploration of the ISM in Galaxies, ed. A. R. Taylor, Tb. L. 
Landecker, & A. G. Willis (San Francisco: ASP), 182 



Lazarian, A. & Pogosyan, D. 2000, ApJ, 537, 720 

Lazarian, A., Pogosyan, D., Vazquez- Semadeni, E., &; Pichardo, B. 

2001, ApJ, 555, 130 
Lazarian, A. &; Pogosyan, D. 2006, ApJ, 652, 1348-1366 
Lazarian, A. & Pogosyan, D., ApJ, in press 

Levy, D., Puppo, G. & Russo, G., 1999, Mathematical Modelhng 

and Numerical Analysis, 33, 547 
Liguori, M., Hansen, F. K., Komatsu, E., Matarrese, S. & Riotto, 

A. 2006, PhRvD, 73, 3505 
Lithwick, Y. & Goldreich, P., 2001, ApJ, 562, 279 
Liu, X.-D. & Osher, S., 1998, Journal of Computational Physics, 

141, 1 

Maccarone, T. J. & Coppi, P. S. 2002, MNRAS, 336, 817 
Mac Low, M.-M. &; Klessen, R. S. 2004, Rev. of Mod. Phys., 76, 
125 

Masahiro, T. &; Bhuvnesh, J. 2004 MNRAS, 348, 897 

McKee, C. F. & Ostriker, E. C. 2007, ARAA, 45, 565-687 

McKee, C. F. &; Tan, J. C. 2002, Nature, 416, 59 

Nordlund, A. K. & Padoan, P. 1999, In Proc. Second Guillermo 

Haro Conf., ed. J. Franco & A. Carraminana (Cambridge Univ. 

Press), 218 

Ostriker, E. C, Gammie, C. F. & Stone, J. M. 1999, ApJ, 513, 259 
Ostriker, E. C, Stone, J. M. & Gammie, C. F. 2001, ApJ, 546, 980 
Padoan, P. & Nordlund, Al999, ApJ, 526, 279 
Padoan, P. & Norlund, A. 2002, ApJ, 576, 870 
Padoan, P., Jimenez, R., Juvela, M., & Nordlund, A. 2004, ApJ, 
604, L49 

Padoan, P., Juvela, M., & Pelkonen, V.-M. 2006, ApJ, 636, LlOl 
Padoan, P. et al. 2007, ApJ, in press 

Passot, T, & Vazquez-Semadeni, E. 2003, A&A, 398, 845 

Peeebles, P. J. E. & Groth, E. J. 1976, A&A, 53, 131 

Scalo, J., Vazquez-Semadeni, E., Chappell, D. & Passot, T. 1998, 

ApJ, 504, 835 
Scoccimarro, R. 2000, ApJ, 544, 597 

Stanimirovic, S., L. St aveley- Smith, J. M. Dickey, R. J. Sault, and 
S. L. Snowden: 1999, The large-scale HI structure of the Small 
Magellanic Cloud. MNRAS302, 417. 
Stanimirovic, S., & Lazarian, A. 2001, ApJ, 551, L53 
Tynan, G. R., Moyer, R. A., Burhn, J. & Holland, C. 2001, Phys. 
Plasmas, 8, 2691 

Vazquez-Semadeni, E., Ballesteros-Paredes, J. & Rodriguez, L. F. 

1997, ApJ, 474, 292 
Vazquez-Semadeni, E. & Garcia, N. 2001, ApJ, 557, 727 
Wiebe, D. S. & Watson, W. D. 2001, ApJ, 549, L115 



